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Exact solutions are found for the Chandrasekhar Page angular equation which results when the 
Dirac equation in a Kerr Newman space time is separated into its radial and angular parts. The 
solutions turn out to be remarkably simple in form while satisfying the asymptotic conditions de- 
duced earlier. The eigenvalues are found to be the square root of the total angular momentum as 
first found by Dirac for flat space; supplemented by a term which is the product of the mass of the 
Dirac particle times the specific angular momentum of the black hole. The additional contribution 
is what is expected from frame dragging. 
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I. INTRODUCTION 



The ordinary differential equations which result from the separation of the Dirac equation in a Kerr (Kerr-Newman) 
space time were first derived by Chandrasekhar and Page The Kerr Newman line element in Boyer Lindquist 

coordinates is given by Q 
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Here M is the black hole mass, Q the black hole's electric charge and a is the hole's angular momentum per unit 
mass. We will only consider the case M 2 > Q 2 + a 2 and work in natural units G = 1, c = 1 and % = 1. As this metric 
is axially symmetric, the Dirac equation can then be separated by writing the four components of the spinor wave 
function in the form, 
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The notation follows Chandrasekhar 's book slightly modified to use uj for the frequency (energy) rather than a. 
The R±i are functions of r alone, while S±i are functions of 0. If this form is substituted into the Dirac equation, 
we get the following angular equations, 
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i5_j_i = (— A + am e cos 9)S_i. 
S_i = (A + am e cos 9)S + i . 



Again for the radial equations it is convenient to use a slightly different variable from that used in 
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we define 



'Department of Physics, Indian Institute of Technology Delhi, New Delhi 110016 India; Electronic address: dilip@physics.iitd.ac.in 



2 



R i = vAi?i to obtain, 
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The separation constant A was shown by Carter and McLenaghan jj] to be the eigenvalue of square root of 
the total angular momentum operator in the asymptotic limit. This was recently confirmed in an analysis of the 
Chandrasekhar separation ansatz by Batic and Schmid 

These two first order equations for the angle J2J|, can be combined into a second order equation, which is known as 
the angular Chandrasekhar Page equation, 
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The equations © and © have been the object of study since their discovery I n the above references, 

recurrence relations for the series solutions as well as series for the eigenvalues have been found. The numerical 
nature of these has meant that the solutions are approximate and there have been disagreements as to the actual 
numerical values [9j- In addition, due to the very large number of parameters; three for the black hole M, Q, a, two 
for the Dirac particle m e , e and two for the dynamical system as a whole u>, m; the possible solution space is very 
large. In the references cited above numerical results have been found only for a few cases. 

The full Dirac equations have also been studied by Finster and others ^3 an d Dafermos who claim that no 
solutions exist for the Dirac equations in this case. For the extremal case M 2 = Q 2 + a 2 , Schmid [l2( has found exact 
solutions. The system has also been studied in integral equation form by Guihao Tian . Many of these workers 
expand the solution in the orthonormal basis formed by the spin weighted spheroidal harmonics |14l Il5l Il6l Il7| | . This 
appears to be the natural basis to expand as in the limit M = 0, Q = 0, the metric eq.Q reduces to the Minkowski 
metric in oblate spheroidal coordinates. The Dirac equation is separable in this coordinate system and the angular 
solutions are spin weighted spheroidal harmonics. If further a = also, these reduce to spin weighted spherical 
harmonics 1 1 lL llq. 

The emphasis in this work is on the requirement of a finite probability density everywhere for the wave function. 
In the next section, for the case m = 1/2, I show that if this condition of finite probability density is imposed, an 
explicit and simple solution is possible only for certain values of the parameters. The eigenvalues are also found in 
this case. 



Based on this solution (|15|) . in the third section the problem is also solved for arbitrary non-zero m. An explicit 
expression for the eigenvalues and the allowed values of the parameters is obtained. The eigenvalue expression is 
very similar to the expression obtained in the case of m = 1/2 but the states are very different in nature. Indeed, the 
solutions, equations l|23|l are so simple in functional form that it is easier to verify them by direct substitution into 
the angle equations © rather than solve the equation. Along with our earlier result 18] for the special case m = 0, 
this completes the solutions of the angular Chandrasekhar Page equation. In the last section we discuss these results. 



II. SOLUTION OF THE EQUATIONS [3] FOR m = 1/2 

For m — 1/2 the angle equations J2Jl become 

S + i = (— A + am e cos9)S_i, (6a) 
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We isolate the part where this equation can possibly be singular by rewriting the equations as 
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The substitutions 
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convert our equations © to 
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These can be easily rewritten as a second order differential equation, 
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The above equation gives a Sturm Liouville problem with regular singular points at x = ±1 and A/(am e ) and hence 
the solutions corresponding to acceptable wave packets can always be expressed in terms of hypergeometric functions. 
We confine ourselves to the case where the eigenvalue |A| > |am e |. The justification for this restriction is given in 
the last section. The power series solution will then converge over the entire allowed range of x, namely [—1, 1]. We 
therefore try a series solution of the form 



n=0 



n=0 



Substituting the above in equation @ the resulting recurrence relations for the series coefficients are, 
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This series does not converge at both the limits ±1 and thus cannot represent a wave packet over all of space. 
Further convergence of the series for T + i alone is not enough as the function of actual interest S + i has a sin (0/2) 
appearing in the denominator, equation iJHJ. To overcome the divergence due to this sine, it is necessary to have 
A 2n — — A2 n +i. Imposition of this condition leads to determination of the eigenvalues A and solutions, which are 
possible only for 
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The solutions for m = 1/2 are 
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These wave packets have their maxima at one pole of the black hole and their minima at the other. 

III. SOLUTION OF THE EQUATIONS |3] FOR ARBITRARY rn = k + 1/2 

For arbitrary m = k + 1/2 the angle equations (JSJ become 
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As in the previous case, we use the substitutions given in equation JJJ and equation JHJ to obtain 
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These can be easily rewritten as a single second order differential equation with regular singular points at x = ±1 
and \/{am e ). Once again we confine ourselves to the case where the eigenvalue |A| > |am e |. The power series solution 
will then converge within the entire allowed range of x, namely [—1,1]. As before, we try a solution of the form, 

On substituting the above in equation (|17fl . the resulting recurrence relations for the series coefficients are, 

nC n — {k + au)C n -\ — (n - 2)C n _ 2 + awC n _ 3 = AD n _i - (flra e + A)D„_ 2 + am e D m - 3j (19a) 
n.D„ + (fc + aw)Z) n _i - (n - 2)Z?„_ 2 - awD n _ 3 = -AC„_i - (am e + A)C*„_ 2 - am e C„_ 3 . (19b) 
Again the second equation of the pair becomes an identity if 



A, 



P{-l) n C n 



P = ±1. 



(20) 



C n then satisfies 



nC n -{k + auu- PA(-l) n ) C„_i + (n - 2 - (am e + AP(-l)™)) C n _ 2 - (aw + am e (-l)"P) C„_ 3 = 0. 



(21) 



This series does not converge at both the limits ±1 and thus cannot represent a wave packet over all of space. 
Further convergence of this series for Ti alone is not enough as the function of actual interest 5 + i has a sin (0/2) 
appearing in the denominator, equation ©. To overcome the divergence due to the sine, it is necessary to have 
= — C^n+i- Imposition of this condition leads to determination of the eigenvalues A and solutions, which are 
possible only for 
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It is interesting to note that unlike the flat space case, only positive values of k are now allowed. That is only 
corotating wave packets are possible. The solutions for m = k + 1/2 are 
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These go over to the asymptotic forms found in |9( at the poles and have their maxima close to the equatorial 
plane. 



IV. DISCUSSION 



The solutions given in equations (|15|l and l|23|) are so simple in form that the it seems surprising that they have 
not been obtained before. The reason appears to be that they are not easily expressed in series composed of the spin 
weighted spheroidal harmonics requiring an infinite number of these terms for their representation. Thus previous 
efforts were using an ill adapted basis and consequently did not obtain an easy convergence. In our case the factoring 
out of the half angle term is responsible for the quick convergence of the remaining part. 

There is a simple physical explanation for the occurrence of the half angles. The action of the spin weight 1/2 
raising and lowering operators [lj, [l5j on these states is as follows. 
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Such states with s > m are explicitly excluded in the definition of the spin weighted harmonics, as they are defined 
as functions of spin and total angular momentum |l4l Il5|. In the language of flat space quantum mechanics, the 
solution expanded in terms of the spheroidal states are given in the basis \J,S). The states as expanded here are 
in the basis \L, S). That is, we conceive of states like i^oo a s a product of a spin and a orbital angular momentum 
state. Both sets of course span the same space and will yield the same results, only the second set happens to be 
more efficient in this situation. 



Consider first the eigenvalues we have found. Our neglect of all eigenvalues |A| < \am e \ is for the following 
reason. In the classical domain, these correspond to trajectories with impact parameter less than a, thus they 
will penetrate deep inside the horizons and are expected to be captured. In quantum terms, the spin con- 
tribution of h/2 alone is many orders of magnitude greater than the even extremal values of a for any particle 
with mass less than a picogram and thus the neglect of states with |A| < |am e | is justified for all known Dirac particles. 

It is interesting to note that the eigenvalues given in equations p4h) and l)22hl are the linear sum of the eigenvalue 
found by Dirac for the flat space central force problem along with a contribution which is the product of the 
mass m e of the particle and a the specific angular momentum of the hole. Due to frame dragging we expect that 
the particle will indeed acquire such an additional contribution to the angular momentum as seen by a distant observer. 

There is one important difference in the eigenvalue spectrum from the flat space case. In that situation 
k = ±(j + 1/2) and for a given j all permitted values of the azimuthal quantum number —j < m < j can occur. 
Here the only value which occurs is m = j. That is the hole ensures that all allowed states with non zero orbital 
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angular momentum are in the equatorial plane and are co-rotating. Further the total angular momentum of the 
Dirac particle is aligned with that of the hole m = j. This is to be expected on classical grounds. 

The separation of successive eigenvalues by two units of angular momentum is also to be expected given the spin 
2 nature of the gravitational field. It is perturbations of this field which will cause transitions between the states of 
different angular momenta. This is analogous to the situation in the case of atoms where the angular momentum 
states are separated by one unit of angular momentum, corresponding to the spin 1 nature of the electromagnetic 
field. This causes the spectrum of states for an atom to contain only states of half integral angular momenta or only 
integral angular momenta. 

Next consider the states. The maxima for all the states with m ^ 1/2 are very close to the equatorial plane 
and converge very rapidly to it with increasing m. Equally interesting is the fact that this displacement from the 
equatorial plane is different for spin up and spin down states. The spin up states have their maxima below the 
equator while the spin down states have their maxima above the equatorial plane. This is to be expected on the 
basis of the magnetic and gravimagnetic interactions of the Dirac particle in the non uniform field of the hole. 

The states with m = 1/2 are exceptional as these alone have their maxima at the poles. Again with spin up at the 
south pole and spin down at the north pole. Along with the confinement of the other states with non-zero orbital 
angular momentum, this is strongly reminiscent of the accretion disk and polar emission of classical particles orbiting 
a black hole. Of course, until the radial Chandrasekhar Page equations are also solved this is only a analogy. 

The radial Chandarsekhar Page equations can be brought to almost the same form as the angle equations studied 
here and thus may be amenable to solution by the same method. 



S. Chandrasekhar, Proc. Roy. Lond. A. 349, 571 (1976). 
D. Page, Phys. Rev. D. 14, 1509 (1976). 

S. Chandrasekhar, The Mathematical theory of Black holes, Oxford University Press, Oxford, UK 1992. 
B.Carter and R.B. McLenaghan Phys. Rev. D . 19, 1093 (1979). 

D. Batic and H. Schmid; e-print gr-qc/0512112 

E. D. Fackerell and CM. Cosgrove; J. Math. Phys., 18 1849 . (1977) 

E. J. Kalnins and W. Miller, Jr.; J. Math. Phys. 33. 286 (1992). 

S.K. Chakrabarti, Proc. Roy. Lond. A. 391, 27 (1984). 

D. Batic, H. Schmid and M. Winklmeier; J. Math. Phys., 46, 012504, (2005), e-print gr-qc/0402047| 

F. Finster, N. Kamran, J. Smoller and S.T. Yau ; Adv. Theor. Math. Phys. 7 . 25 (2003) e-print |gr-q c/9905047 
M. Dafermos; e-print gr-qc/0210025 Commun. Math. Phys. 238 (2003), 411427. 

H. Schmid; e-print gr-qc/0207039 
Guihua Tian; e-print gr-qc/0411087 

E. T. Newman and R. Penrose; J. Math. Phys. 7, 863 (1966). 

Goldberg. J. N., A.J. MacFarlane, E.T Newman, F. Rohrlich, E.C.G. Sudershan; J. Math. Phys. 8, 2155, (1967). 

A. Staruszkiewiez; J. Math. Phys. 8, 2221, ( 1968) 

E. Berti, V. Cardoso and M. Casals, ; e-print gr-qc/0511111 
D. Ranganathan; e-print gr-qc/0306090 

P.A.M. Dirac, The principles of Quantum Mechanics, Oxford University Press, Oxford, 1991. 



